Two loop results from the derivative expansion of the blocked action 
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A derivative expansion of the Wegner-Houghton equation is derived for a scalar theory. The 
corresponding full non-perturbative renormalization group equations for the potential and the wave- 
function renormalization function are presented. We also show that the two loop perturbative 
anomalous dimension for the 0{N) theory can be obtained by means of a polynomial truncation in 
the field dependence in our equations. 
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In field theory the description of physical observables 
at energy scales that can be several orders of magnitude 
below the cut-ofF scale, is one of the most relevant issues. 
This problem is studied by means of the renormahza- 
tion group (RG) and, in particular, the Kadanoff- Wilson 
formulation of the RG M , based on the blocking proce- 
dure, is probably the most preferable starting point for 
any application. The Wegner-Houghton equation 0] is 
a specific realization of the blocking procedure that re- 
lates the action defined at the scale k to the action at 
the scale k — Sk through a functional integration over the 
infinitesimal momentum shell 5k, for (5fc — > 0. This differ- 
ential equation is exact, as shown in B since it contains 
the whole effect due to the integration of the fast modes. 
However it has no known analytical solution and one has 
to resort to approximations in order to deal with a more 
tractable problem. 

The most promising approach consists in a derivative 
expansion where the action is expressed as a sum of terms 
containing an increasing number of field derivatives. The 
lowest order corresponds to a momentum independent 
projection of the Wegner-Hougton equation which is then 
reduced to a renormalization group equation for the local 
potential H. To the next order, 0{d^) terms are included 
and one gets two coupled differential equations for the 
potential V and the wave-function renormalization func- 
tion Z. In yj and, more recently, in g, these equations 
are studied numerically to investigate fixed points and 
determine critical exponents. Some applications of the 
0(9^) truncation are considered in |^^. 

Here we focus on the connection between perturba- 
tive calculations and the various approximations of the 
Wegner-Houghton equation. The exactness of the full 
equation implies that one should be able to recover the 
perturbative results order by order in the coupling con- 
stant at least in some approximation scheme. For in- 
stance if one considers a polynomial expansion in the 
field 4> for the potential F of a scalar theory and re- 
duces the equation for ^ to a set of differential equa- 
tions for the various polynomial coefficients (m^,A, ...), 
it is very easy to check that the equation for the quartic 
coupling A yields the perturbative one loop /3-function, 
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provided that one neglects all couplings corresponding to 
non-renormalizable powers of (p (0®, 0^, ...) and takes the 
ultraviolet (UV) limit where the mass is negligible com- 
pared to the running scale pOl. However, going beyond 
one loop is no longer straightforward. 

In this paper we specifically consider the anomalous 
dimension 77 for a 0{N) symmetric scalar theory, which 
gets the first non-vanishing contribution at two loop, and 
show how the O(A^) term is obtainable from the Wegner- 
Houghton equation. In fact the two loop anomalous di- 
mension and /3- function have been derived in [[ll] p^ , 
by making use of the average action formalism; however 
this result is not entirely satisfactory since the authors, 
by employing that formalism, are forced to introduce a 
smooth cut-ofF, namely an exponentially decreasing func- 
tion, in order to integrate out the UV modes. Instead, in 
the following it is shown that the very simple and intu- 
itive scheme which yields the one loop /^-function can be 
extended to deduce r/ at two loop level, just solving the 
differential equation obtained for a particular truncation 
of the polynomial expansion of Z, not resorting to any 
smooth cut-ofF. 

Moreover, starting from the two loop anomalous di- 
mension, we compare expansions of the action having 
different nature and leading to different equations for Z. 
We show that the plain derivative expansion introduced 
in [|l^,|5| has the correct UV behavior, whereas the one 
suggested by Aitchison and Fraser llfl] is misleading in 
the UV region. As a consequence we can discriminate 
between the two approaches and choose the first one in 
order to consider any further application. 

Let us consider the Euclidean iV-component scalar the- 
ory (a, 6 — 1,2, ..TV and (j)^ = 4'i'P^) in Four dimensions, 
and, according to our choice of considering a derivative 
expansion of the action, truncated at 0{d^), the action 
at the scale A reads 

SM ^ J d^x(^Z''''d,,cbad^<l>b + v) (1) 

where the potential V and the wave-function renormal- 
ization function Zab are polynomial in (f>i. In the follow- 
ing the 0{N) symmetry of the theory is not required, 



except when explicitly indicated. The blocked action 
5'fe[$], which is the action at the scale fc < A and which 
is a functional of the field $(x) that, by construction has 
Fourier components only up to the scale fc, is obtained 
by functional integration from 5a [(/i]- Up to terms that 
are not relevant for the Wegner-Houghton equation, the 
functional integration yields 



Skm^SAm + ^Tv'lniK,, 



S^JKo) 



(2) 



where Kij = {S'^SA/64>iS4>j)\^=q, and Tr' is the trace of 
the internal degrees of freedom and of the modes within 
the shell fc — A, and finally SijKo is the free propagator 
with mass toq- 

The Wegner-Houghton equation is the differential 
equation for Sk [^] obtained by differentiating Eq. (||) with 
respect to fc, with the boundary condition fixed by S'a [</']■ 
By double-differentiating Eq.(Q) we get 
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where the notations Ai^j\ — [Aij -\- Aji)/2 and Aij^k = 
dAij/d^k have been used. A comment is in order: in de- 
riving Eq. (0) no integration by parts has been performed 
and thus no surface term in the action has been neglected; 
indeed these terms, which are irrelevant in the equation 
of motion, supply a higher order (in the weak coupling 
perturbative expansion) contribution to Eq.(g). 

By employing the conventions S^{x ~ y) = {x\y), 
{d/dXfj,)6^{x-y) = -i{x\p^'\y) = -{d/dy^,)5'^{x-y) and 
the commutation rules [p'",/(a;)] = id^'f{x), [p'^,f{x)] — 
2id'^f{x)p^ — d'^f{x), Kij{x,y) can be cast in the oper- 
atorial form {x\Kij\y) with 



+ \Zab,^,d^<^^d''<^'' - d^ (%b), ,9^$'' 



(4) 



where we have ordered all coordinate operators to the 
left of the momentum operators and 
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(5) 



Next step concerns the evaluation of the trace in Eq.(0). 
An expansion of the trace in derivatives of the fields 
is required in order to determine an evolution equation 
for V and Zab^ and such a suitable expansion has been 



already derived in ||lJ,|T5[. It is obtained by defining 
Kij{u) = Kij + u6ij and Ko{u) = Kq + u, where u is 
a fictitious "mass" , then deriving with respect to u and 
finally expanding the derivative of the logarithm. By in- 
troducing 



r(M) = 2Ti"'ln {Kij{u) - 5ijK(3(u)) , 
the required term in Eq.(y), r(0), is given by 



(6) 



r(o) 



du—T(u) 
du 

duTr'{K-/{u)-S,,K-\u)) (7) 



Note that since the trace is performed on a finite momen- 
tum shell no IR or UV divergence affects Eq. (M). Then 
we split Kij{u) = Aij + Bij -t- C^ with 

Aij = ^(jj)P^ + V;.y B,j = Z(ij) {p^ - p^) 
Cy = iXi" ^d^'<^bP^. 

A, = \Zab,^,d^^''d^'<^'' - d^ (%,), ,9'^$'') (8) 

where the eigenvalue p^ of the momentum ( p^\p) — 
p^\p)) has been introduced. Finally K~Mu) in Eq. ffl) 
can be expanded according to 



K-^{u) = A-^ -A-^{B + C + D)A~^ 
+A-^{B + C + D)A-\B + C + D)A-'^ 



(9) 



For a detailed analysis of the validity of this operatorial 
expansion we address to the Report of R. Ball ||l^ . 

In order to deduce the evolution equations for V and 
Zab, it is sufficient to retain respectively all terms with 
no or two derivatives. Actually C and D contain respec- 
tively one and two field derivatives; then in the trace in 
Eq. (0) it is convenient to order all operators with all B's 
on the right side since, by definition, B\p) = 0. Therefore 
B appears in the trace only inside commutators and, due 
to the commutation rules, each commutator containing 
B provides one or two derivatives. It follows that the 
truncation displayed in Eq. (|9|) is sufficient for our pur- 
pose. As shown in p5|jl6[ the trace is then reduced to 
an integration over the spatial coordinates as well as over 
the momentum modes in the k — A shell. At this level it 
is allowed to integrate by parts and neglect the surface 
terms since derivatives of 6 are no longer present; making 
use of (a = l/(167r2) ) 



d^p 



(2^)4 
we finally get 

Sk = Sa + 



PtiP" 9{P^ 
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(10) 
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where the trace only concerns the matrix indices i and 
j, J indicates that the integral is restricted to the shell 
P - A? and 

+A-^M''A-~^R^A-^ + A-i(i?" + S"' - X'')A-H'I''A-^ 
+-A-^^S''A~^X^A-^p^ + -A-^X^A-^ZA-^M^A-^p^ 

-A-^S-'A-^ZA-^M^A-^p^ - -A-^X''A-^X''A-^p^ 

+A-HrA-^ZA-^ZA-Hl''A-^p^ 

--A-^M^A-^ZA-^X^A-^p^ 
2 ^ 



(12) 



V^ki") and Q''\i = Z''\m and 



with Affe,'^ = (^(fcO'V 

Rk'^i = Zffe") , and Ski" = Zkf". It is understood that 

in the r.h.s. of Eq. (p^), the omitted lower indices are 
contracted, according to the given sequence, in such a 
way that, with a and b fixed, each term of the sum, as 
well as T'^^, is a matrix with lower indices i and j. 

The equation for the potential is then recovered from 
the piece not containing field derivatives whereas the 
d^^ad'^^b coefficient gives the equation for Zab- 

For a general A'^-component theory it is not possible to 
commute the matrices in Eq. (11) and this makes the 



integration of the variable m, which appears only in A~^, 
a hard task. Beside a direct numerical approach there are 
special cases where Eg. (pi]) can be drastically simplified. 
For instance in ref. [|l5| it is shown how to perform the 
integral through a change of variables if Zab — Sat in the 
r.h.s. of Eq. (|ll|); however this is equivalent to neglect 
all Zab contributions to the scale evolution of Zab which 
would be governed by the potential only. 

In the following we consider two special cases in which 
Eq. ( [ll| ) can be easily treated. The first one consists in 
checking that for the 0{N) symmetric theory, the pertur- 
bative two loop anomalous dimension is recovered from 
Eq. (lO). To this aim it is sufficient to truncate the 
general expansion of Zab in powers of $ up to 0($i$j) 



Zab = ZoSab + Zi^'^Sab + Z2^a^b 



(13) 



and the potential to order 0($^): V = m'^^^ /2 + X^'^/A\. 
In order to calculate the anomalous dimension, the func- 
tion ZQ{k) must be determined. Therefore one has to re- 
tain in Eq. (|ll| ) the numerical coefficient of Sabd^^'^df^^'' 
only, which corresponds to set $ = in the coefficient 
of Sab9'^^°'di^^^. In this case the w-integration is easily 
performed. According to Eqs. (pi]), (O) and by assum- 
ing, as boundary conditions of our problem, Zab — 6ab at 
the scale A, we get 



zq 
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(14) 



In order to solve Eq. dl4) the corresponding equations 
for zi , Z2 are needed. However as it will be clear at the 



end, instead of considering the full equations for these 
two couplings, it is sufficient to retain the contributions 
proportional to the quartic coupling squared A^. This is 
achieved by taking Zab = ^ab in the r.h.s. of Eq. (|l] 



and collecting terms proportional to 5ab^ d^^"^d^^^ and 
$a$65^$°a^$''. We get 
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Z2 
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X / dp^ p^ 



3(^2+^2)3 4(p2+TO2)4 



(15) 



Eqs. (pH), ([15| ) can be solved perturbatively order by 
order. We start with all couplings fixed by the boundary 
conditions at A, i.e. za — 1, zi = Z2 = 0, put these values 
in the r.h.s. of the equations and solve them, obtaining 
the first order solution for zq, zi, Z2. Then the first order 
couplings can be put in the r.h.s. of the equations to get 
the second order solution and so on. By this procedure 
we get to the first order, in the UV region (fc^ >> m^): 
zo = 1, zi =: (aA2)/(54 k"^), za = (iV + 6)(aA2)/(108 P). 
Effects 0(fc^/A^) are neglected. To the second order for 
Zq we get zo = 1 - {a^X'^iN + 2)/18)ln(fc/A) and the 
corresponding anomalous dimension {t = ln(fc/A)) 

,._!,„ « = ;!!i!(i^ + 0,A», ,16, 

is in agreement with the standard result ||l^. It is easy 
to realize now that, since the neglected terms in the dif- 
ferential equations for zi and Z2 are proportional to non- 
renormalizable couplings, which are zero to the lowest 
order, they can only provide O(A^) contributions to zq 
and to the anomalous dimension. 

This result not only shows that the particular trun- 
cation adopted here still retains the features of the two 
loop perturbative calculations but mainly shows that the 
presence of the sharp cut-off k in Eq. ( pi] ) does not af- 
fect quantities such as the anomalous dimension which 
is related to a two loop diagram characterised by over- 
lapping divergences. It follows that all difficulties raised 
in |11|] due to the use of a sharp cut-off and the neces- 
sity of introducing a smooth cut-off in the framework of 
the average action, disappear when the sharp cut-off is 
embedded in the formalism developed above. 

As a second example we consider Eq. (^l]) for the 
A^ = 1 theory. In this case there are no more matrices and 
the u-integration is trivial. Thus the Wegner-Houghton 
equation in this case reads (here A = Zk"^ + V" and the 
prime indicates derivative with respect to $) 






A 



(17) 



dZ^ 
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„ ,4/Z" Z'A 
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3^3 3A3 4A4 

These two non-perturbative coupled equations are suit- 
able for a numerical investigation of the fixed point struc- 
ture of the theory as explained in ^Q . 

We shall now discuss the difference between Eqs.(p^, 
(|l8|), and the analogous equations which are derived by 
resorting to the logarithm expansion proposed in [16|. 
This is an expansion of different kind, as stressed by [17|, 
based on the presence of "small" fluctuations of the field 
<&(x) around a constant background $o; this allows to ex- 
pand the logarithm in Eg . (B ) up to some chosen power 
of the fiuctuation field [|l6|,|7| . According to this expan- 
sion the differential equations for V and Z have been 
derived in Q. Although near the Gaussian fixed point 
the linearised form of Eq. dia) , and of Eq. (7) of ref. M are 
the same, some differences appear at the quadratic level. 
To illustrate this point in a specific case it is sufficient to 
examine the equations for Z when only the effects arising 
from the potential are taken into account. In this case 
one gets ( §1) 



, dZk 
' dk 



-ak\V'y 



V" 



(fc2 + V"Y 



while from Eq. ( [l8|) one instead finds 



dZk 
' dk 



ak\V"'Y P + AV" 
6 (fc2 + V"Y 



(19) 



(20) 



One notes that the UV behavior is different in the two 
cases. It is remarkable that this difference is responsible 
for the non-zero anomalous dimension obtained above. 
Conversely the expansion due to Aitchison and Eraser 
leads to a zero anomalous dimension to O(A^) and the in- 
clusion of higher order terms in the derivative expansion, 
such as 0{d^) or 0{d^), does not change this result. One 
can also check that the r.h.s. of the two equations behave 
very differently near the region of the mass gap, thus pre- 
dicting different scaling behavior even at the crossover. 
Only when one integrates the modes, in the independent 
mode approximation, from A to A: << m, Eqs. (19^ (120) 
differ by 0{k'^ /m^) and as already noted in P, |15| , |17| , |l9| , 
by integrating up to fc = 0, one gets in both cases 



Zk=o — 1 



{V"'f 
1927rV' 



(21) 



On the basis of our result in Eq.(pii[), we argue that the 
derivative expansion illustrated above is more suitable 
for deriving a quasi-local approximation of the Wegner- 
Houghton including 0(9^") terms in the blocked action. 

In conclusion we have derived two non-perturbative 
coupled differential equations for the potential V and the 
wave-function renormalization function Z , by making a 
derivative expansion of the blocked action and employ- 
ing a sharp cut-off. These equations when compared to 
the usual perturbative result show a good UV behavior 
providing the correct two loop anomalous dimension for 
the 0{N) symmetric scalar theory. 
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